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ABSTRACT 

We have detected the Sunyaev-Zel'dovich (SZ) increment at 850 /im in two galaxy clus- 
ters (CI 0016 + 16 and MS 1054.4 - 0321) using SCUBA (Sub-millimetre Common User 
Bolometer Array) on the James Clerk Maxwell Telescope. Fits to the isothermal (3 model 
yield a central Compton y parameter of (2.2 ± 0.7) x 10~ 4 and a central 850 /im flux of 
A7o = 2.2 ± 0.7 mJy beam -1 in CI 0016. This can be combined with decrement measure- 
ments to infer y = (2.38±g;ff) x 10~ 4 and w poc = 400±i|°°km s" 1 . In MS 1054 we find a 
peak 850 /mi flux of A/ = 2.0±1.0mJybeam _1 and y = (2.0 ± 1.0) x 10 -4 . To be success- 
ful such measurements require large chop throws and non-standard data analysis techniques. 
In particular, the 450 /im data are used to remove atmospheric variations in the 850 /im data. 
An explicit annular model is fit to the SCUBA difference data in order to extract the radial 
profile, and separately fit to the model differences to minimize the effect of correlations in- 
duced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can 
be used to measure the SZ increment in massive, compact galaxy clusters. 

Key words: cosmic microwave background - cosmology: observations - galaxies: clusters: 
general - large-scale structure of universe - methods: observational - submillimetre 



1 INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect can be used as a powerful 
probe of cosmology. Because the intensity of the SZ distortion is 
virtually independent of redshift, SZ measurements provide an av- 
enue to study clusters and their peculiar velocities to any distance. 
Additionally, the SZ effect allows independent determination of 
various cosmological parameters. For example, because the SZ ef- 
fect's amplitude is proportional to the electron density along the 
line of sight through the cluster, n c , while the X-ray amplitude is 
proportional to n 2 ., SZ effect data can yield distances to clusters in- 
dependently from the cosmological distance ladder. This informa- 
tion can then be used to p rovide estimates of the Hubble constant 
H fe.g. lReese et all2000h . 

Th e physics of the SZ effe ct is qui te sim- 
ple ISunvaev & Zel'dovich! Il97l IBirkinshawl 1 19991 
ICarlstrom. Holder & Reesd l2002h . Cosmic microwave back- 
ground (CMB) photons can scatter off hot (T ~ 10 7 K) electrons 
in a plasma. A statistical net gain of energy of the photons from the 
electrons occurs, producing a characteristic distortion in the CMB 
as seen through the plasma. This distortion produces a decrement 
in the CMB's temperature below about 200 GHz, but an increment 
in the CMB's temperature above this frequency. 

The most important use of SZ increment measurements lies in 
the determination of the spectral shape of the SZ effect in clusters 
of galaxies. This shape is the sum of the thermal and kinetic effects 



together with other sources of emission within the cluster. The ki- 
netic SZ effect is a change in the apparent CMB temperature as seen 
through the cluster due to the peculiar motion of the cluster relative 
to the CMB. The kinetic SZ effect may be present in any given clus- 
ter and may manifest itself as either a temperature increase or de- 
crease, although it never dominates the thermal effect for expected 
cluster velocities jBirkinshawlll999l) . Knowledge of the cluster's 
properties obtained at several wavelengths can allow separation of 
these two effects. If the kinetic effect can be isolated, the cluster's 
peculiar velocity along the line of sight can be determined. Statis- 
tics about large scale flows at a variety of redshifts provide strong 
constraints on the dynamics of structure formation. In principle, de- 
tailed measurement of the SZ spectral shape can also constrain the 
cluster gas temperature through the relativistic correction ( Rephaeli 
1995, Birkinshaw 1999, Colafrancesco, Marchegiani, & Palladino 



High frequency measurements can also help separate contam- 
inants such as dusty galaxies which may exist in millimetric mea- 
surements, since the dominant sub-mm point sources differ from 
those seen at longer wavelengths. 

Measurement of the Sunyaev-Zel'dovich (SZ) effect incre- 
ment in galaxy clusters has been an elusive and important goal 
of observational cosmology for over a decade. While measure- 
ments of the SZ effect decrement a re common IBir kinshawl 
Il999l ICarlstrom. Holder & Reesell2002l) . measurement and study 
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of the increment has been less successful. This is largely be- 
cause current instruments generally lack the sensitivity to mea- 
sure small amplitude, extended, positive signals in the sky. Also, 
other positi ve flux sources in a cluster field can c onfuse SZ mea- 
surements lLoeb&RefregIeilll99l iBlaml Il998t) . However, ex- 
periments work ing near the null of the th ermal SZ effect like 
SuZIE jHolzapfel etall 1 9971) . and ACBAR freterson et a?]|2002t) 
are beginning to yield results. The balloon-borne experiment 
PRONAOS has had limited success in SZ increment measure- 
men ts iLamarre et alll998l) . and ot her detections are als o claimed 
(e.g. lKomatsu^tTL^999Tr Recentlv. |Benson et alJl2003l) have per- 
formed this type of measurement in six galaxy clusters, and find 
that measurement of any given cluster's peculiar velocity is diffi- 
cult. 

Measurement of the SZ increment is difficult for a variety of 
reasons, including the amplitude of the atmospheric signal at sub- 
mm wavelengths. 450 /im and 850 /im observations take advan- 
tage of windows in the opacity of the sky at these wavelengths, 
but successful observation of any astronomical source in the sub- 
mm regime still requires non-trivial atmospheric subtraction tech- 
niques. Additionally, because the SZ effect is a small amplitude 
signal extended across a relatively large region of the sky, sensi- 
tive instruments which can sample a range of spatial frequencies 
are required. Care must be taken when reconstructing the spatial 
shape of the SZ effect from data recorded as differences between 
two or more telescope pointings. As point sources are expected to 
be a major contaminant in the sub-mm regime, high resolution is 
also required to adequately understand their effect in a given field. 
One of the instruments currently most capable of performing obser- 
vations of the SZ increment is the Sub-millimetre Common User 
Bolometer Array (SCUBA; iHolland etailll999h . This instrument 
provides enough sensitivity and angular resolution to detect and 
remove point sources (including lensed background sources) from 
the data. SCUBA also samples large enough regions of the sky to 
constrain the shape of the SZ increment in moderate to high redshift 
clusters. However, even with these characteristics, measurement of 
the SZ increment with SCUBA is not trivial. 



2 OBSERVATIONS 

The bulk of the data discussed here was obtained on November 10 
and 11, 1999 at the James Clerk Maxwell Telescope (JCMT) at the 
summit of Mauna Kea, Hawaii. The JCMT combines a high, dry 
site with a 15 metre dish and sensitive instruments. SCUBA, at- 
tached on the Nasmyth platform of the JCMT, allows simultaneous 
observations with 91 450 /im-band bolometers and 37 850 /im- 
band bolometers. Because the 450 /J,m band observes very little SZ 
increment flux, while the expected peak of the SZ signal occurs 
near 850 /J,m, the 450 /im data can act as a monitor of atmospheric 
emission in this experiment. 

SCUBA is sparsely sampled spatially, which is undesirable 
in this experiment. 'Jiggling' to fully sample the array consists of 
small (~ 6 arcsec) changes in the position of the beam in a regular 
pattern. The detectors integrate in each jiggle position for 1.024 s. 
Although 64 positions are required to fully sample both arrays, this 
experiment used 16 positions because it was only necessary to fully 
sample the 850 fim array. Moreover, executing the jiggle pattern as 
quickly as possible is more conservative in terms of residual sky 
fluctuations. 

Variations in the signal produced by atmospheric (sky) noise 
and thermal offsets from the instrument are removed by rapid chop- 
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Figure 1. Our chopping strategy in relation to the SZ increment. The solid 
curve is an isothermal (i profile (discussed in Section 3.3) with parameters 
9 C = 50 arcsec and /3 = 0.75. The profile's peak has been normalized 
such that AIo = 1.0. The dashed lines show the SCUBA array size at the 
'source' position marked 'A', and the two reference positions, 'B' and 'C 
The filled areas marked 'd' show the FWHM of one bolometer, and indicate 
the bolometer's position over an observation. SCUBA jiggle mapping gives 
data in the form h = c^a — (ciB + d-c)/2 for each bolometer. Because it 
is sensitive only to differences in SZ flux between the source and two refer- 
ence beams, h contains much of the information about the SZ effect ampli- 
tude, although this is the same information contained in the shape within the 
array size. However, the JCMT does not have high enough resolution to rely 
on the differenced SZ shape alone. Using the 850 fim array average (that 
is, the mean of every bolometer's h) to remove the atmospheric emission at 
each time step subtracts a flux about the same size as h, which destroys the 
information about the SZ effect amplitude in this data. 

ping together with slower telescope nodding (e.g. I Archibald et all 
2002). Rapidly chopping the secondary mirror faster than the rate 
at which the sky is varying removes much of this atmospheric sig- 
nal. Nodding the telescope's primary mirror to a second reference 
position removes slow sky noise gradients and instrument noise. 
These motions produce a chop pattern such that the two reference 
beams are each observed half as often as the source beam (Fig. 1). 
In the standard SCUBA mapping mode, the chop occurs at about 8 
Hz and the nod occurs at about 1/16 Hz. Although a chop throw 
in the range 50 to 100 arcsec is standard, this experiment requires 
a chop that does not put significant SZ flux in the reference beams 
(Fig. 1). For this reason, an extremely large JCMT secondary mir- 
ror chop size of 180 arcsec was used. 

89 ks of data were taken, split between several fields. The CSO 
225 GHz r-meter was not operational for either observing shift, 
so we relied on skydips to measure the atmospheric optical depth, 
T850, in our data set. Both nights had relatively constant tsso, but 
the first evening had a lower average (tsso — 0.15) than the second 
(rsso — 0.3). Uranus was used as a calibration source throughout. 

Three clusters and three blank fields were observed during 
this run (see Table 1 for details), referred to as CI 0016, MS 0451, 
MS 1054, Blank 02, Blank 09 and Blank 25 hereafter. 23 ks of ad- 
ditional data were collected with an aluminum reflector masking 
the aperture of SCUBA. These data are used as null tests to study 
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SCUBA's behaviour under both optical and electrical influences 
(blank sky data) and electrical influences only (reflector data). 

Additional data were taken at the JCMT on October 2, 3, and 
6, 2002 in the same manner as the 1999 run. Both CI 0016 and an- 
other blank field, Blank 21, were mapped. In addition to the jiggle 
mapping, photometry of possible point sources in the CI 0016 field 
were performed. The photometry used a 9 point jiggle pattern with 
1 arcsec spacing and a 60 arcsec chop throw. Bolometer G9 of the 
SCUBA array was not operating properly during this run and its 
data are not used here. A total of 13 ks of data were taken in 2002; 
how they are combined with the earlier data is discussed below. 



3 GENERAL ANALYSIS 
3.1 Data Preparation 

The data are deswitched, flat-fielded and corrected for atmo- 
spheric extinction with SURF (the SCUBA User Reduction Facil- 
ity, Jenness & Lightfoot 1998). Extinction correction is performed 
by calculating the airmass at which each bolometer measurement 
was made, and then multiplying by the zenith sky extinction at the 
time of the measurement to give the extinction optical depth along 
the line of sight. Each data point is then multiplied by the expo- 
nential of the optical depth to give the value that would have been 
measured in the absence of the atmosphere. The photometric cali- 
bration of deep SCUBA fields like those presented here is known 
to be accurate to about 15 per cent teorvdl200l lArchibald etafl 
2002), although ~ 2/3 of this uncertainty is due to error in the 
absolute calibration of the planets' flux. 

It is well documented that the deswitched 450 /im 
and 850 /im time streams are highly correlated, and that 
this correlation stems primarily from atmospheric noise 
(e.g. IJenness. Lightfoot. & Holland Il998t iBorvs et alJ Il999l 
Archibald et al. 2002). The standard method of removing at- 
mospheric noise from an array is to subtract the array average 
from each bolometer at each time step. However, most of the 
information about the SZ flux in this measurement is contained 
in the difference between the source and reference beam fluxes 
(see Fig. 1) rather than in the shape of the SZ emission across the 
field of view. Therefore, removing the 850 /im array average at 
each time step would drastically reduce the inferred SZ emission. 
However, because the SZ effect is small at 450 /im, using that 
band's data should mitigate this problem. A least-squares linear 
fit of the 450 /im array average to the 850 /im array average at 
each time step is used to obtain a linear coefficient Ci and an 
offset Co (hereafter referred to as the (450 /im) method). These fit 
parameters and the 450 /im array averages are then used in place 
of the 850 /im array average to subtract the atmospheric noise at 
each time: 

S? 50 (t) = Rl 50 (t) - (d (R 450 (t)) + Co) , (1) 

where Sf 50 (£) is the corrected data at time step t for bolometer 
6, Rf 50 (t) is the raw, double-differenced and extinction corrected 
data at time step t for bolometer b, the angled brackets denote an 
average over the whole array, and the superscripts denote wave- 
length. The fits are performed separately for each (approximately 
20 minute) observation. C\ removes the atmospheric emission on 
scales of about 1 s, while Co removes the atmospheric average over 
periods long compared to the telescope nod time. The value of Ci 
and Co do not vary by more than 1 per cent over an observation; 
this is consistent with the results of B orvs et alJ ll9991. C\ is typi- 
cally 0.6 and Co is of order 1 /iV in our data set. The linear Pearson 



correlation coefficients of the 450 /im and the 850 /im time streams 
are on average 0.95 and never less than 0.90 for these observations. 
The average ratio of the variance of S$ in the time stream of the 
(450 /im) atmospheric subtraction method to the variance of the 
850 /tm time stream in the standard SURF atmospheric subtraction 
method is 1.15. This implies that the (450 /tm) method is slightly 
inferior for removing atmospheric noise. However, this approach 
is required for the analysis of SZ data, since it allows retention of 
information on angular scales up to that of the chop throw. 

Of the seven target fields, three are unusable for this exper- 
iment in some way. A fundamental requirement of this analysis 
method is high enough signal to noise ratio that it is possible to 
determine the location of point sources in a field and remove them. 
However, the Blank 02 field had only a short integration time (see 
Table 1), hence it is not included in the results of this experi- 
ment. Blank 22 was a poor choice for a control field, because 
it contains larg e point sources which corrupt it for our purposes; 
Chapman et al. 1 2001) examine this field in detail. The MS 0451 
field contains a sub-mm bright gravitationally lensed arc serendip- 
itously discovered when com bining this data with previous obser- 
vations (Chapman et al. 2002). Although an interesting target in its 
own right, the ~ 12 mJy bea m" resolved arc dominates the sub- 
mm emission in this cluster iBorvs et alJl2003l) . making measure- 
ment of the SZ increment impossible with these data. 

3.2 Point Source Removal 

Recognition and removal of point sources, especially unbiased han- 
dling of lensed background sources, is important in this experiment 
jLoeb & RefregieJl99l.lBlairll998l) . We adopt the following pro- 
cedure. The measured flux differences are binned according to the 
spatial position of the 'source' beam. The standard method of pre- 
senting a SCUBA image is to convolve the binned data set with the 
telescope point spread function (PSF). Maps made by this method 
are shown in Fig. 2 in order to indicate the locations of possible 
p oint sources and to ex hibit their consistency with the same fields 
in lChapman et alJ l2002t) . However, these images are not used as 
part of the quantitative analysis. 

Instead, the binned flux differences are convolved with a 
wavelet similar to a 'Mexican Hat' function: 

w>=* ( i -4)™ p (^)- (2) 

Here 9 is an angular distance and <j w is a parameter charactering 
the wavelet's width. The normalization tpo is set such that the flux 
from a point source placed in a blank field is returned identically at 
8 — 0. This wavelet is more useful than the telescope PSF in this 
analysis because it returns an amplitude which is insensitive to the 
average value in the map, which is not expected to be zero in an SZ 
experiment. The wavelet's width was chosen to be <r w = 6.3 arc- 
sec, slightly smaller than the a corresponding to the FWHM of the 
telescope, because this value maximizes the signal to noise ratio in 
simulated measurements. After convolution with this function, the 
positive and negative extrema in the maps are the locations of pos- 
sible point sources in the source and reference beams respectively 
(Table 2). 

The philosophy of this experiment's point source removal 
scheme is not to develop a list of point sources per se, but rather 
to remove the flux contributed by any point-like source. In the sim- 
ple method utilized here, possible point sources in a field are re- 
moved by adopting a single flux level above which sources are 
found. The limiting flux level is estimated based on atmospheric 
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Table 1. Experimental target field parameters, including field names and coordinates, on source integration times, and dates of observation. 

Target Field Type RA (2000) DEC (2000) Integration Observation 

Time (ks) Dates 

Blank 02 Blank Field 2:22:59.9 +4:59:57 2.5 11 Nov. 1999 

Blank 09 Blank Field 9:00:00.0 +10:02:00 9.4 10,1 1 Nov. 1999 

Blank 21 Blank Field 21:30:00.0 +16:30:00 7.2 03,06 Oct. 2002 

Blank 22 BlankField 22:17:25.1 +0:12:59 8.8 11 Nov. 1999 

CI 0016+16 Cluster Field 0:18:33.7 +16:26:04 12.2 1 1 Nov. 1999, 02,03 Oct. 2002 

MS0451-03 Cluster Field 4:54:11.5 -3:00:52 15.1 10,1 1 Nov. 1999 

MS 1054-03 Cluster Field 10:56:57.4 -3:37:24 15.1 10,1 1 Nov. 1999 



CI 0016+16 



MS 1054-03 




00 h 18 m 40 s 




10 h 57 m 04 s 



Blank 9h field 



10° 00' 30 




09 h 00 m 04 s 



08 h 59 m 56 s 



16° 29' 00 




21 1 30" 06 



Figure 2. Maps of our target fields using a standard map-making procedure (similar to the SURF map-making algorithm). These maps are not used at any point 
in our analysis, except to serve as a check that large point sources are not present in the data. Contours show the point source cutoff levels in each field (quoted 
in the text). The la rms values of these maps are: 2.9 mjy in CI 0016; 3.3 mJy in MS 1054; 3.6 mjy in Blank 09; and 3.9 mJy in Blank 21. The fluxes and 
(possible) locations of point sources found and removed from later analysis are given in Table 2. 
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Table 2. Candidate point sources which are removed. 



Field 


Flux (mjy) 


RA (2000) 


DEC (2000) 


CI 0016 


10 


0:18:32.7 


+16:25:18 


MS 1054 


-11 


10:56:52.5 


-3:37:27 


MS 1054 


-12 


10:56:53.0 


-3:37:45 


Blank 09 


13 


9:00:02.4 


+10:01:02 



extinction, integration time and the noise equivalent flux density. 
For each field, we use the following thresholds: 9 mJy in CI 0016, 
9 mJy in MS 1054; 11.5 mjy in Blank 09; and 13 mjy in Blank 
21. It is found that t he point source list matches that found in 
IChapman et alj J2002) well. The estimated variance given by this 
method and the actual mean variance of the pixels agree to 10 per 
cent or better in each field. 

By determining the reference beam positions at which most 
of the negative flux originates, it is found that the 'two' negative 
point sources in MS 1054 are, in fact, probably caused by one 
bright source. The source appears to be two sources because the 
MS 1054 field was observed only at the beginning and the end of 
the Nov. 1 1 1999 shift. The field rotated during the interval, thereby 
extending the single source into 'two' apparent sources. Our best 
estimate of this source's true position is 10:57:03.5 ±2.0 seconds 
in rig ht ascension and —3:39:00 ±60 arc seconds in declination. 
Chap man et alj 120021) find another source of similar amplitude in 
this cluster which did not happen to fall in any of our beams. The 
presence of these sources is consistent with the source count mod- 
els (e.g. lSmail et all2002l) . 

The effect of greater than 3a positive sources is removed 
directly from the time series using the measured amplitude, the 
telescope pointing history and the telescope's PSF. A negative 3a 
source is equally likely to be in either of the two reference beam 
positions, so the same approach is not available. Hence, data as- 
sociated with detected negative sources are simply excised from 
further analysis. 

Monte Carlo simulations have been performed to develop an 
understanding of the effects of point sources on this data. These 
simulations are f ully explained in Zemcov, Newbury & Halpern 
(2003) (hereafter ZNH), and the application of the results found 
there are discussed in Sections 3.3 and 5. 

3.3 Model Fitting 

Because the SZ effect is extended on the sky, the simplest way to 
isolate it in a cluster is to fit the data to a model. Since the signal 
to noise ratio in any given pixel is poor for this data set, a simple 
isothermal /3 model is adequate: 

/ fl 2\ (l-3ffl/2 

AI D (0) = AI? (l+^J • O) 

Here, AIq is a measure of the magnitude of the increment at the 
centre of the cluster, 9 is the angular distance from the centre, and 
8o and P are parameters characterizing the cluster. These parame- 
ters have already been determined for these clusters from combina- 
tions of X-ray and SZ decrement data (see Table 3). Although there 
is some uncertainty in the values of j3 and 8q, the error this causes 
in the inferred Al® is negligible at this experiment's level of preci- 
sion. We thus fix the angular dependence and focus on determining 
the amplitude alone with our SCUBA data. 

Because the SZ increment is non-zero at 450 /im and the 



Table 3. Isothermal (i profile fit parameters. 



Cluster 


do 


P 


z 


Reference 


CI 0016 


42.3" 


0.749 


0.55 


Reese et al. (2000) 


MS 1054 a 


66.0" 


0.75 


0.83 


Jeltemaetal. (200 n 


Blank 09, 21 


50.0" 


0.75 




Generic Values 



a Jeltema et al. (2001) discuss the difficulty of applying an 
isothermal (3 profile to MS 1054, which is a double-cored 
X-ray cluster. 



Table 4. SZ model fit results. 



Target 


AT" J 3 (mjybeam- 1 ) 


cr(AI^) (mjybeam- 1 ) 


CI 0016 


1.5 


0.5 


MS 1054 


1.5 


0.7 


Blank 09 


-0.4 


0.7 


Blank 21 


1.1 


1.0 



450 /im array average is used to remove residual sky noise, the su- 
perscript 'D' is used as a reminder that the value quoted is actually 
a difference between the 850 /im signal and the 450 /im array aver- 
age. Retrieving the 850 /im signal (denoted AJq 50 ) from AJq 3 will 
be discussed further below. 

Directly fitting the data to an isothermal /? model differenced 
according to the telescope pointing history is the most straightfor- 
ward way to determine the peak SZ increment flux in a given field. 
Both the A/g 3 and the la error in Aig found from fitting to the 
differenced isothermal j3 model are listed in Table 4. The error in 
AJq is found via the correlation matrix of the linear fit. a(AI^) 
is the square root of the second diagonal element of this correlation 
matrix. The fit correlation matrix has very small off diagonal ele- 
ments, meaning that the fit is not correlated with the (meaningless) 
voltage offset value. 

To test the hypothesis that much of the information about the 
SZ intensity is contained in the difference Sf 50 (t), its average is 
calculated over the array at each time step. The mean of these val- 
ues should be greater than and be detected with approximately the 
same or slightly less significance as the detection given by the fit to 
the isothermal /3 model. It is found that the means of the CI 0016 
and MS 1054 fields are approximately 1.5<r and 2a above zero, 
respectively. We conclude that merely using the shape of the SZ 
distortion across the SCUBA field of view underestimates its true 
amplitude. 

Note that intensities, AJsze, are generally quoted in mjy per 
JCMT beamsize, which is related to temperature change via the 
derivative of the blackbody distribution: 

where x — hv/kTcMB is the dimensionless frequency, Tcmb = 
2.725 K, ATszb is the change in CMB temperature due to the SZ 
effect, and SI = 5.74 x 10 -9 steradians per beamsize. 

It is also desirable to make a crude 'map' to check the 
radial profile of the data. A matrix inversion technique uti- 
lizes all of the available data to make a m ap, including in- 
forma t ion about the off-source pointings ( e.g. iTeemark & Bunn 
I l995t IWrieht Hinshaw & Bennett! Il996l IStompor et all l200l 
Hinsh av^eUdl2003l) . This method involves dividing the sky into a 
vector of pixels, creating matrices based on pointing and flux data, 
and the inversion of a single matrix to create a map. The matrix 
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inve rsion itse lf is performed using a Singular Value Decomposition 
(e.g. lPress et alll992t) . 

Here, the pixel model of the sky is chosen to be the bright- 
nesses of annular rings centred on the cluster, as this conforms to 
the approximately circular symmetry expected in the SZ signal. The 
rings have 40 arcsec widths to track the radial variation with rea- 
sonable signal to noise ratio. Once the fluxes of the annuli including 
their errors and correlations have been estimated, they can be fit to 
the isothermal f3 model. The fit is performed using both data and 
model with a weighted mean of 0. This ensures that the model's DC 
term does not falsely increase the measured increment. A general- 
ized x 2 statistic is determined, where x m is a pixel's isothermal [3 
model value, and Xd is a pixel's data value. 

X = (x m -x d )i Cy 1 (x m - X d )j. (5) 

Our observational strategy produces non-negligible correlations be- 
tween the annular pixels which are given in Cy, the pixel-pixel 
covariance matrix. Explicitly, dj = (P T N _1 P) , where P is 
our pointing matrix and the noise is N. N is assumed to be diag- 
onal because removing the fit to the 450 /im array average effec- 
tively acts as a pre-whitening filter. These correlations need to be 
carefully included when the SZ amplitude is determined using this 
annular pixel model and an isothermal j3 profile. Fig. 3 shows the 
annular pixel fluxes as well as the best fit profile (Equation (3)) us- 
ing the parameters listed in Tables 3 and 4. The two methods yield 
very similar results, although we find that the direct fitting method 
is more sensitive. 

Generally, the figure of merit for the goodness of fit is that x 2 
should be approximately the number of degrees of freedom in the 
data. However, source confusion may have an effect on the appli- 
cability of this figure of merit. To test how well x 2 estimates the 
goodness of fit, Monte Carlo simulations similar to those presented 
in IZNHl are performed in order to determine the probability that the 
cluster field fits are consistent with the absence of an SZ increment. 
The direct fit to the isothermal /3 model differences method is used 
to find AJq 1 in simulations containing realistic point sources and 
with noise levels set to be the same as those in this experiment, but 
with no SZ increment in the fields. Point sources are removed at 
the same level as in this experiment; Fig. 4 shows the results of 200 
such simulations. It is found that the fits associated with the high- 
est 10 per cent of the x 2 statistics always yield drastically incorrect 
increment values because of point source contamination. We there- 
fore exclude the worst 10 per cent of the fits. Fig. 4 also shows the 
SCUBA results for both clusters, each of which pass the x 2 cutoff. 
The results overlap, as the fits for these fields give almost the same 
increment value. If the SZ increment were absent from CI 0016 and 
MS 1054, the probability of obtaining the AJ found is less than 
approximately 0.01 for either cluster. 

The x 2 values found for CI 0016 are 3.7 for 6 degrees of free- 
dom via the seven pixel fit method, and a x 2 of 178402 for 178406 
degrees of freedom using the direct fit to model differences. We 
conclude that the fits are fully consistent with the presence of an 
increment in this cluster. 

Our confidence in the parameters found for MS 1054 is less 
than that for CI 0016. Based on its x 2 of 13.4 for 6 degrees of 
freedom from fitting the annular pixels to an isothermal /3 model, 
this cluster would normally be flagged as a poor fit. The fit may 
be poor because the isothermal (5 model does not describe the dis- 
tributi on of MS 105 4's intracluster electron gas particularly well 
Ijeltema et all200ll) . However, the fit may become tolerable if the 
nature of the double-cored electron gas distribution is taken into ac- 
count. Nevertheless, the reduced x 2 from fitting to the isothermal f3 




Figure 4. Results of simulations to determine the effect of confusion on 
these measurements. This plot shows the cumulative probability of finding 
an increment value A7o if no increment were present, P(AI > Ai"o), as a 
function of fit amplitude ATq found from fitting to the differenced isother- 
mal j3 model. The simulations use realistic noise including point sources 
and no SZ effect. Point source removal has been implemented (see Section 
3.2 for details). The values determined for the two clusters in this study are 
shown by the dashed vertical line. This function is not symmetric about zero 
because gravitational lensing preferentially brightens positive point sources, 
thereby biasing the point source removal scheme. 



model differences is 258438 for 258204 degrees of freedom in this 
cluster. The reduced x' 2 statistic associated with this fit falls well 
below the worst 10 per cent of reduced x 2 statistics in the Monte 
Carlo simulations discussed above. The fit is therefore acceptable 
based on this criterion. 

The parameters determined for the blank fields are also given 
in Table 4. From the fits to annular pixels, the Blank 09 field has a 
X 2 of 11.7 for 6 degrees of freedom and the Blank 21 field has a x 2 
of 13.6 for 6 degrees of freedom. The x 2 from fitting to the model 
differences are poor in both fields. Both fields also seem to exhibit 
unremoved, confused point sources. The simulations o f lZNHl show 
that fields with no SZ increment form a well-defined locus on the 
Io~X 2 plane, and both of the blank fields in this experiment would 
be inside the confused regime and hence be rejected as SZ incre- 
ment candidates based on their x 2 ■ 

For the purposes of comparison, the line A/o = OmJy 
beam -1 is fit to the data using the direct fit method. The Ax 2 (= 
X 2 — Xmin) values for these fits are 6.7 for CI 0016, 6.6 for 
MS 1054, 0.5 for Blank 09 and 2.0 for Blank 21. These values are 
in good agreement with the significance of the detections. 

Although the MS 0451 and Blank 22 fields were found to be 
unsatisfactory for this experiment earlier in our analysis, this analy- 
sis pipeline also rejects these fields. Both fields have AIo consistent 
with mjy per beam, and poor x 2 due to the number and location 
of sources within them. This is an objective check that retrieving 
an SZ increment from data heavily polluted by point sources is dif- 
ficult. 
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Figure 3. Data points determined via matrix inversion, together with best fit isothermal (3 profiles found from a direct fit for the 4 target fields shown as ±1<t 
confidence regions by the grey band. Note that the isothermal /3 profiles are not best fits to the data points shown but rather use the model in Table 3 fit directly 
to the SCUBA differences. The annular points are used to assess the fields' radial behaviour and to determine the effect of point sources in the beams. As 
matrix inversion yields data with a mean of zero, the mean of the best fit model is added to the data in these plots. The zero level is shown as a dashed line. 



4 INDIVIDUAL CLUSTER ANALYSIS & RESULTS 

In this section our data are used to determine the 850 /im SZ in- 
crement values in CI 0016 and MS 1054. A maximum likelihood 
method is used to determine the best fit Compton y parameter 
(which parametrizes the integrated pressure in the cluster along the 
line of sight) and, in CI 0016, the peculiar velocity of the cluster 
along the line of sight. 



4.1 CI 0016+16 AI^ 50 Determination 

IWorrall & Birkinshawl <2003l) find that CI 0016 has T c ~ 
9keV, implying that it has a relativistic electron population 
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which may give a significant correction near the positive 
peak of the SZ spectrum . Ther efore, the corrections found by 
lltoh. Kohvama & Nozawal 1 19981) . which are good to fifth order in 
k^T e /m e (? , are used to determine the correct shape of the thermal 
SZ distortion. These corrections, which represent a 10-20 per cent 
change to the flux at 850 /im for standard cluster parameters, are 
applied to all of the spectral calculations, and are implicit in the use 
of the term 'thermal' SZ effect below. 

Because the SCUBA filters have a finite bandwidth, we need 
to determine the effective central frequencies of this measurement. 
This calculation must be performed because the SZ spectrum dif- 
fers from the grey-body spectra SCUBA usually observes. To do 
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Table 5. Equivalent u, A for SZ spectra. 



Filter 


V 


c/v 




850 fim 


347 GHz 


865 


fim 


450 /im 


652 GHz 


460 


fim 



this, the 'equivalent frequency' is defined as: 

_ _ / °° vF(y)S(v)dv 



/~ F{u)S{u)dv 



(6) 



where F(v) is the actual SCUBA filter's response and S(v) is the 
SZ effect spectrum over the region of interest. The results of this 
calculation are listed in Table 5; we conclude that the frequency 
shift from the nominal values are negligibly small. 

Foreground dust emission could also contaminate the mea- 
surement of the SZ increment in these clusters. ISO 1 data are used 
to estimate the dust contribution. The telescope pointings are used 
to sample the 180 and 90 ISO map of the cluster region in 
exactly the same way as the SCUBA data. Errors for each mea- 
surement are determined from the error map and incorporated into 
the simulated time stream. These data are then fit to the model dif- 
ferences. The x 2 values imply that the fits are very poor; this is 
expected because any dust emission in the region (Galactic cirrus, 
extragalactic point sources, or the cluster itself) is unlikely to be 
distributed like the isothermal (3 model. However, the dust may still 
contribute to the SZ signal and so must be removed. The dust con- 
tribution is determined by fitting a modified black body of the form: 



Id(v, T d ) 



.3+7 



ghi//fc B T d 



(7) 



Here, the parameters are the spectral index 7, dust temperature Td, 
and normalization a, which is determined using both of the ISO 
data points weighted by their errors. 7 is taken to be 2.0 and Td 
to be 20 K, parameters w hich are typical for 7 50-detected cirrus 
at these wavelengths (e.g. lLamarre et alll998l) . Using this model, 
the dust emission at 450 is determined to be Id (660 GHz) = 
0.07 ± 0.04 mjy beam -1 for CI 0016. This value is not very sen- 
sitive to variations in the parameters 7 and Td. The fit is corrected 
for dust emission by the addition of this value, which is a few per 
cent correction. 

Because our analysis scheme subtracts the 450 fim flux, a 
method of retrieving A7q 50 from A7° is required. The actual 
quantity subtracted from each 850 fim bolometer's time stream is 
the average of the 450 array at each time step. This average 
contains the SZ signal sampled as described in Fig. 1, with A7q 50 
being the peak isothermal (3 model amplitude. The average value of 
the 450 fim array is determined by sampling the model SZ profile 
in the same way as in this experiment. A map is made based on 
this data, and the average of this map in the field of view is con- 
structed. It is found that there is a linear relationship between the 
input SZ amplitude and the average of the differenced map, given 
by A7f v ° = AAI$ 50 , where A is a linear coefficient and A7l 5 v ° is 
the average of the difference map over the field of view. This means 



that A7q 50 can be found via: 
A7 n 8 "" 



A7 I5 + J 4A7 450 , 



(8) 



1 We have used archival data from the Infrared Space Observatory (ISO), a 
European Space Agency mission with the participation of ISAS and NASA 
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Figure 5. Sources of extended emission in the CI 0016 field. The dashed 
line shows the corrected thermal SZ effect determined from both our 
850 fim point and low frequency decrement measurements. The inset shows 
the low frequency data. The dot-dash line shows the dust emission de- 
termined from 750 measurements (asterisks). The best fit kinetic effect is 
shown as a dotted line. The sum of the thermal and kinetic effects is shown 
as a solid line. For reference, the position of 450 fim (dotted vertical line) 
and zero flux (horizontal solid line) are shown. Our 850 fim point is shown 
as a diamond, and has been corrected for both the 450 fim subtraction and 
dust emission. 



where A depends only on the shape of the isothermal /3 model. 
A — 0.72 in CI 0016, including the gain difference between the 
850 fim and 450 fim arrays. 

To determine the best fit Compton y parameter and pecu- 
liar velocity u pcc , we utilize a likelihood function which searches 
the y - v pcc parameter space for its maximum. Previous decre- 
ment measurements at 15 GHz I Gr ainge e t al] |20 02f). 2 GHz 
(Hugh es & Birkinshawj[l993) and 30 GHz jReese et allEoOch are 
used alone to determine a best-fit y value, and also in conjunc- 
tion with our measurement to help constrain both the thermal 
and kinetic SZ effects in CI 0016. Alone, these values yield y = 
(2.38 ± 0.12) x 10~ 4 in this cluster, assuming v pcc = km s" 1 . 
Our value alone yields y = (2.2 ± 0.7) x 10~ 4 under the same 
assumption. This measurement is therefore consistent with earlier 
results. 

We find that the 850 fim SZ increment value for CI 0016 is 
A7,f = 2.2 ± 0.7 mjy beam" 1 . This corresponds to 0.38 ± 0.12 
MJy sr" 1 or ATcmb = 1.2 ± 0.4 mK in thermodynamic temper- 
ature units. Fig. 5 shows the corrected 850 fim data point, along 
with the radio and 750 results. Also plotted are the best fit thermal 
SZ effect and dust emission spectra, and the position of the 450 fim 
band-centre for reference. 

In addition to constraining y, the amplitude of the kinetic SZ 
effect is also constrained by this likelihood analysis. Using both our 
measurement and the three decrement measurements, it is found 
that the marginalized best fit parameters are y = (2.38±o;34) X 



European Space Agen 
iKessler et alll99M) . 



10 and v-r, 



400± 



1900 km s" 1 in CI 0016. Assuming v pe 



0, these points yield y = (2.38 ± 0.11) x 10" . Fig. 5 also shows 
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the inferred kinetic SZ effect from the 850 ^im data. This result 
is therefore only a weak constraint on the line of sight peculiar 
velocity. 

A kinetic SZ effect in CI 0016 may be mimicked by a primary 
fluctuation in the CMB, because the spectral distributions of the 
two effects are the same. However, the primary CMB anisotropies 
at the scales of this ex periment, £ > 50 00, are about lfiK rms 
Ihu & Dode lson 2002), iBorvs et aDll999l) . This corresponds to a 
±0.005 mjy beam -1 signal, which is negligible at our level of pre- 
cision. 

4.2 MS 1054 03 A/ ( f Determination 

The value of y for MS 1054 is found by determining the differ- 
ence J 850 - Al 450 for the SZ distortions defined by 1.0 x 10 -5 < 
j ^ 1.0 x 10 -3 , where A = 0.81 for this cluster. The y val- 
ues giving differences most closely matching the data and la error 
points are the best fitting y parameter and error bars. This method 
yields y — (2.0 ± 1.0) x 10~ 4 for this cluster, including rela- 
tivistic corrections. These y values correspond to an increment of 
A/, 850 = 2.0 ± 1.0 mjy beam" 1 , or AT C mb = 1.0 ±0.5 mK. 

Because little other relevant data exist for this cluster, fur- 
ther analysis is not possible. The dust contamination cannot be 
estimated becau se there are no ISO data, and the available dust 
maps fcchleeel. Finkbeiner & Davis 1998), which are of poor an- 
gular resolution, provide little spectral information. However, it ap- 
pears that the dust contamination is at a similar level to that found 
in CI 0016, and hence is probably negligible. Although one pub- 
lished estimate of the thermal decrement magnitude already exists 
for MS 1054 Ijov et al 1200 II) . a worthwhile estimate of its kinetic 
effect cannot be made because no explicit y parameter or ATcmb 
is available. 



5 DISCUSSION 

In principle, measurements of the SZ increment, in combination 
with lower frequency measurements, can yield a great deal of infor- 
mation. For example, the line of sight peculiar velocity and possibly 
even the cluster temperature can be determined via measurement of 
the SZ increment. However, measurement of the SZ increment with 
SCUBA is difficult. There are three essential components which 
combined to make this experiment a success. 

The first component is a carefully planned observational strat- 
egy designed to maximize SCUBA's sensitivity to an SZ signal. 
Using a large chop throw to reduce the tendency to subtract SZ sig- 
nal and chopping in azimuth so that sky rotation reduces the effect 
of point source contamination are both parts of this. Also, using the 
450 fim array average to remove the spurious sky signal, thereby 
keeping SZ information on scales larger than the array size, max- 
imizes the available signal in the data. Controls such as the blank 
fields, and to a lesser extent the reflector data, are required to check 
for systematic effects in both the instrument and analysis methods. 

The second necessary component of this experiment is the use 
of custom-designed software to handle the data analysis, which is 
particularly important for fitting the differenced data directly to a 
model. 

The third important component of this analysis is the removal 
of point sources in order to reduce their contaminating effect. Be- 
cause sub-mm bright, high z background sources are relatively 
common in cluster fields, observations of the SZ increment could 
suffer prohibitively from the effects of point sources. The basic 



principle we employ is that point sources affect only a few pixels 
and hence can be distinguished from an extended SZ effect. Our 
data have a confusion limit of approximately 0.8 mjy beam -1 for 
SZ signals. Integration to SCUBA's fundamental confusion limit of 
~ 0.5 mjy beam -1 rms (which would require weeks of integra- 
tion time) gives an SZ confusion of ~ 0.5 mjy beam -1 . Our sim- 
ulations show negligible reduction in the confusion limit of these 
observations whe n point sources with amplitude less than about 5 
mjy are removed l ZNH). This limit would require approximately 
2 shifts of integration time. Unfortunately, decreasing the confu- 
sion limit floor below 0.5 mjy beam -1 would require higher an- 
gular resolution observations. A corollary to this is that SZ incre- 
ment detections from other experiments with lower resolution than 
SCUBA, such as SuZIE and ACBAR, may be confusion dominated 
if unsupported by higher resolution data. 

These measurements confirm the SZ effect in 2 galaxy clusters 
at an amplitude consistent with the measured decrements. A weak 
constraint on the kineti c effect in CI 0016 i s found which is con- 
sistent with the result in lBenson et alJ J2003). We quote a marginal 
2a SZ increment detection in MS 1054, but cannot do much more 
with the data at this point. 

It is clear that large reductions in the uncertainty of the value 
for the kinetic effect velocity in CI 0016 will not be achievable be- 
cause of source confusion with this type of instrument. This means 
that SCUBA is not capable of determining the amplitude of the 
kinetic effect in any specific cluste r, given the e xpected level of 
SZ-derived peculiar velocities ("e.g. ISheth & Diaferidl200ll) . It is 
likely that the only way to use SCUBA for constraining peculiar 
velocities is through a statistical survey of many clusters. 

Nevertheless, we have demonstrated that, with sufficient care, 
it is possible to measure the 850 /im increment at the JCMT. Such 
SCUBA observations of the SZ effect at 850 fim can be combined 
with measurements at other frequencies to study individual clusters 
in more detail. Follow-up observations with high resolution sub- 
mm instruments are likely to become an important part of forth- 
coming SZ cluster surveys. 
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